Spontaneous plaquette formation in the SU(4) Spin-Orbital ladder 
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The low-energy properties of the SU(4) spin-orbital model on a two-leg ladder are studied by a 
variety of analytical and numerical techniques. Like in the case of SU(2) models, there is a singiet- 
multiplet gap in the spectrum, but the ground-state is two-fold degenerate. An interpretation in 
terms of SU(4)-singlet plaquettes is proposed. The implications for general two-dimensional lattices 
are outlined. 

PACS numbers: 75.10.Jm, ll.30.-j, 75.40.Mg 

The properties of Mott insulators with orbital degen- 
eracy is attracting a lot of attention with the increasing 
evidence that this degeneracy can have many other conse- 
quences apart from the standard cooperative Jahn- Teller 
effect. One of the possibilities that seems to be realized in 
LiNiC>2 JO is that the additional orbital degree of freedom 
prevents the system from ordering in both the orbital and 
spin channels. It was suggested in a recent paper by Li 
et al. [Q that this might occur if spin and orbital de- 
grees of freedom play a very symmetric role, like in the 
SU(4) symmetric version of the Kugel'-Khomskii model 
p] defined by the Hamiltonian: 



As we shall see, exact diagonalizations suggest that the 
ground state is a two-fold degenerate plaquette solid with 
gapped multiplet excitations. The important step for- 
ward though is that analytical results can be obtained in 
both the weak and strong rung limits finally putting this 
plaquette picture on very firm grounds. 

The SU(4) spin-orbital model on a ladder is defined by 
the Hamiltonian 



H = J|| ^ l!Q (2s*i,Q-s*i+i,a + |)(2Ti,Q.r i+ i jQ + |) 
+ J± Ei(2si,i.Si,2 + |)(2f i ,i.r ii2 + |), 



(2) 



H = ^2 Jij(28i.Sj 



(i) 



Such a Hamiltonian is indeed a good starting point for 
LiNiC>2 due to the local symmetry and the strong Hund's 
rule coupling, but its properties are only beginning to 
be understood. The fundamental difference with SU(2) 
models stems from the fact that it takes at least 4 sites to 
make an SU(4) singlet. For the ID version of the model, 
which is fairly well understood both at zero and 
finite temperature |8|,|9| , this shows up as a four-site peri- 
odicity of the correlation function. In 2D lattices, it was 
argued by Li et al. ||Jl^] that the system might prefer to 
form local SU(4) singlet plaquettes in the ground state 
rather than developing long-range order. While exact 
diagonalizations (ED) of the model on a square lattice 
indeed support this conjecture the lack of analytical 
results in any limit prevents one from drawing definite 
conclusions. 



where a site on the two-leg ladder is described by its rung 
number i and its chain index a — 1,2, Si_ a is a spin one- 
half operator at site (i,a) and Ti jCt is an isospin one- half 
corresponding to the orbital degree of freedom on the 
same site (see Fig. pi). 
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FIG. 1. The spin-orbital ladder. 

In this Letter, we have adopted another strategy and 
decided to study the simplest lattice in which plaque- 
ttes might form, namely the two-leg ladder (Fig. 1). 
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FIG. 2. Scaling of the singlet-multiplet gap (+) and 
the singlet- singlet gap (X). The dashed line is a fit with 
A sm = AJ^j + Be~ N ^. It is clearly consistent with a non-zero 
value of A sm in the N — ► oo limit, whereas the behavior of 
A ss is consistent with a zero value in the same limit. 

Isotropic ladder - We start by considering the 
case that is closer to 2D, namely the isotropic limit 
J± = J|| — J. Taking advantage of all symmetries (trans- 
lation, rung-parity cind SU(4) quantum numbers 5^ ^., 7~j^ 
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and sr t z ot = J2i s i T i)> we have obtained the low-energy 
spectrum on clusters with 8, 12 and 16 sites with Lanc- 
zos ED using periodic boundary conditions in the chain 
direction. The results can be summarized as follows. For 
all clusters, the ground state is an SU(4) singlet, and the 
first multiplet excitation is at relatively high energy. Be- 
sides, a plot of this energy gap A sm as a function of 1/N 
(see Fig. 2) strongly suggests that this gap remains in 
the thermodynamic limit. In addition to this multiplet 
excitation there is always one low-lying singlet inside this 
gap (2 in the special case of N=8). The splitting between 
this excited singlet and the ground state A ss is plotted 
in Fig. H as a function of 1/N. Although it is difficult 
to draw definite conclusions with only 3 sizes, these re- 
sults strongly suggest that this gap vanishes, and that 
the ground state is two-fold degenerate in the thermody- 
namic limit. 

Fig. U shows the dispersion of the low-lying states for 
16 sites. Two important facts are to be noticed here. 
First, the ground state and the next singlet lie in the 
k = and k = n sectors respectively. Second, the disper- 
sion has a local minimum at k = 7r/2, which announces 
the soft mode at k = tt/2 found in the chain limit where 

j± = o H|. 
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FIG. 3. Dispersion of the low-lying states in the N = 16 
site ladder. The two lowest-lying states with k — and k = n 
are SU(4) singlets. There is a local minimum at k = tt/2. 



All these results can be qualitatively interpretated in 
terms of plaquette coverings of the ladder [p| jic| , pd|l . The 
plaquette is the ground state of the four-site SU(4) spin- 
orbital system. It is the smallest SU(4) singlet one can 
build with a system of degrees of freedom in the funda- 
mental (d — 4) representation of SU(4). It undergoes 
extremely strong fluctuations which minimize the energy 
by link to E n /N — — J. It is thus a very stable object and 
can be used to describe the physics of many realizations 
of this model. For a larger system with N = 4p, p G N , 
one can build tensor-product states as coverings by such 
plaquettes in which the system is 'tetramerized'. 
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FIG. 4. Plaquette coverings of the 5(7(4) ladder. 



In the case of the two-chain ladder, the number of 
such plaquette coverings is two, as shown in Fig. ^a. 
These two coverings differ by a translation by one lat- 
tice spacing along the direction of the ladder. They are 
rung-symmetric, thus having a + rung-parity. A sym- 
metric and an antisymmetric linear combination of these 
states can be built, giving one k = and one k = n, 
+ rung-parity SU(4) singlet state. These quantum num- 
bers agree with our ED results. In the special case of the 
N = 8 ladder, which has the same topology as a cube, 
there is an extra covering, corresponding to the third 
pair of cube faces that can be occupied by each of the 
two plaquettes (see Fig. ||b). So the plaquette picture 
predicts 3 low-lying states in this special case, again in 
agreement with our ED results. Besides, if the ground 
state is a product of singlet plaquettes, multiplet exci- 
tations require the breaking of a plaquette, with a finite 
energy cost equal to 2 J minus a correction due to the 
derealization of this defect, again in agreement with our 
ED results. 

So all the basic features of our numerical results are 
qualitatively reproduced by this simple plaquette picture. 
The following strong coupling approach provides more el- 
ements to support this tetramerization picture. 

Strong COUPLING - We now turn to the strong rung 
limit J± 3> J||. When J|| = 0, the ground state is ob- 
tained by putting each rung in one of its 6 ground states. 
The ground state of a rung can be thought of as the 
6-dimensional irreducible representation of SU(4), or as 
the set of states (spin singlet x orbital triplet) and (spin 
triplet x orbital singlet). To perform strong coupling 
analysis, we thus have to determine the effective Hamil- 
tonian that will lift the degeneracy in this 6 Arrung -fold 
degenerate subspace. 

To first order in Jii , we need only to consider the cou- 
pling between two adjacent rungs. Denoting by (12) and 
(34) the sites of two adjacent rungs, we can actually cou- 
ple them in two equivalent ways: 1 to 3 and 2 to 4 (Hi) 
or 1 to 4 and 2 to 3 {Hz). To first order, the effective 
Hamiltonians corresponding to Hi and H% can be for- 
mally written: 



rreff 

Hl.2 



i,3 



)v-i 2 {j\ 



(3) 



where the sum over i,j runs over the 36 states of the 
Jii = limit. Now, to go from Hf to iff : we ) us ^ have 
to exchange sites 3 and 4. But this transforms any ket 
\i) (respectively bra in the sum of Eq. (j|) into — \i) 



2 



(respectively since this permutation just changes 

the sign of the singlet and leaves the triplet invariant. 
Given the form of the effective Hamiltonian, we thus have 
fjctf H c ^ 

Now the sum of these Hamiltonians Ho = Hi + Hi is a 
very simple operator because each site is coupled to both 
sites of the opposite rung. In terms of the 15-dimcnsional 
vector A of each rung, whose components are the gener- 
ators of SU(4), it can be written 



H Q = -j-[A 12 .A u ] + J,, 



(4) 



with A12 = A\ + A2 and A34 = A3 + A4. As shown in 
Rcf. pH , this Hamiltonian can be rewritten in terms of 
Casimir operators as: 



Ho = 4J,|Ci 



234 



•4J||(d 



J\\ 



(5) 



So the spectrum obtained when coupling two irreducible 
representations of dimension 6 consists of three levels 
with degeneracy 1, 15, 20 and with energy — 4(Jii + J±), 
— 4J|| and — 4(J|| — J±) respectively. The spectrum of 



Ho is thus linear in J\ 



II- 

Ho 



So Ho S — H , and since 
= Ho-, we reach the conclu- 



Hf = Hf and H 1 
sion that Hf = \H Q . 

A more pedestrian way to reach this conclusion con- 
sists in calculating the spectrum of H\. For small J\\ , the 
levels are linear in Ju, and we have checked numerically 
that the splittings and degeneracies of Hi correspond to 
Ho I '2 when Ju is small. Back to the ladder, the first-order 
Hamiltonian thus writes up to a constant 



H 



eff 



J± 



ij,tot 



(6) 



where A^tot = An + A42. Ho is thus nothing but the 
SU(4) Hamiltonian for the A i tot rung degree of freedom. 
In other words, the effective Hamiltonian is the ID SU(4) 
model in the antisymmetric 6-dimensional representation 
(see Fig. ||). This situtation is analogous to going from 
the 5=1/2 SU(2) spin ladder with ferromagnetic rungs 
to the S = 1 SU(2) spin chain. 
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FIG. 5. (a) The strong-coupling low-energy effective hamil- 
tonian for the SU (4) spin-orbital ladder. ( x ) stand for d — 4 
IR degrees of freedom while (*) stands for d = 6 IR degrees 
of freedom. 

This simple form of the effective Hamiltonian has very 
interesting consequences. The d = 6 representation of 
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SU(4) is a self-conjugate and antisymmetric representa- 
tion of an SU(N) group with N even. It thus falls in 
the cases where the Lieb-Schulz-Mattis- Affleck theorem 
]l2| , ^3[ states that the SU(4) Hamiltonian should either 
have a non-degenerate ground-state followed by gapless 
excitations or have a degenerate ground-state. Affleck, 
Arovas, Marston and Rabson |Q have shown that the 
ground-state is a two-fold degenerate singlet, and breaks 
translation invariance. More precisely, the two ground- 
states are spontaneously dimerized, nearest-neighboring 
sites forming SI/ (4) singlets either between neighbors 
(2n, 2n+l) or between neighbors (2n + 1, 2n + 2). Above 
these ground states there is a gap to magnon-like or 
soliton-like excitations [|15[|14| 1. DMRG calculations 0] 
confirmed this picture of dimer order with short-range 
spin-spin correlations (A-spin correlation length of the 
order of the lattice spacing). 

This strong-coupling regime is very similar to the 
physics we have characterized numerically in the interme- 
diate coupling regime and strongly supports our plaque- 
tte interpretation. First of all, the spectrum has the same 
properties in both cases: The two plaquette-states break 
translation symmetry, have a short correlation length, 
and the first excitation has to be built breaking one pla- 
quette, thus leading to a gap. Besides, a tetramerization 
of a ladder is equivalent to a dimerization in terms of 
rungs. We now come to the weak-coupling regime to 
show how the situation sets in when two SU(4) gapless 
chains are coupled to form a ladder. 

Weak coupling - The weak coupling approach pro- 
ceeds in an analogous way as in the SU(2) ladder. In 
the absence of interchain coupling ( J± = 0) the Hamilto- 
nian (1) describes two decoupled SU(4) spin chains and 
is exactly solvable by the Bethe ansatz [Q. The sys- 
tem is gapless and the low energy physics is described 
by six (three for each chain) massless bosons and is con- 
trolled by the fixed point Hamiltonians of two decoupled 
Wess-Zumino-Novikov-Witten (WZNW) SU(4) X models 
H with central charge c = 3 + 3 = 6. As in the SU(2) lad- 
der, the strategy to tackle with the weak coupling regime 
is to look at the stability of the infrared fixed point with 
respect to the interchain coupling. To this end one needs 
the low energy expressions for the SU(4) spin densities 
in terms of the WZNW fields which has been obtained in 
Refs. iJlTl 



3aR + 3aL 



c/2a j^-A 



H.c. 



-1) 



(7) 



where are the 15 SU(4) spin densities of chain in- 
dex a = 1,2 with components J^rl at As = 0, at 
2kp = 7r/2ao, and at 4A;^ = tt/clo- The uniform part 
of the spin density is the SU(4) spin current with scaling 
dimension do = 1. The other oscillating parts, and 
are WZNW primary fields with scaling dimensions 



3 



d-2k F = 3/4 and d&k P = 1 and transform respectively into 
the fundamental and the antisymmetric two-rank tenso- 
rial representations of SU(4). With these results, one can 
obtain the low energy effective Hamiltonian of the weakly 
coupled SU(4) spin ladder 



H 



2ttv 



2 



eff 



= 1 



32B,3\L + n l n 2) 



(8) 



where we have dropped as usual the marginally irrele- 
vant current-current in-chain interactions as well as the 
interaction between the current of the two chains with 
the same chirality that renormalizes the spin velocity. 
The interacting part of Eq. (||) has two contributions. 
One comes from the uniform and 4k f parts of the spin 
densities (0). It is marginal with scaling dimension 2. 
The other contribution, which stems from the 2kF spin 
densities, is a strongly relevant perturbation with scaling 
dimension 3/2 and thus governs the low energy behavior 
of the model. As an immediate consequence we conclude 
that a gap A ~ Jj_ opens in the spectrum. The deli- 
cate point however is whether or not some gapless modes 
survives in the infrared. This issue can be investigated 
by means of the Abelian bosonization of the SU(4) spin 
densities (0). Using the results of Ref. jlT] we have ex- 
pressed the effective Hamiltonian (j^) in terms of the six 
bosonic fields that describe the ultraviolet fixed point. 
The resulting bosonized Hamiltonian is too lengthly to 
be reproduced here but it can be shown that all degrees 
of freedom are massive. 

At this point, the physically relevant question is 
whether the "plaquette" picture drawn from the strong 
coupling analysis survives at weak coupling, and in par- 
ticular whether the nature of the low lying excitations 
at strong coupling changes as the interchain coupling is 
reduced. In the SU(2) ladder, the nature of the low en- 
ergy spectrum is the same in both limits and is captured 
by the weak coupling approach. In this case it has been 
shown in Ref. |l8| that the effective Hamiltonian sepa- 
rates into two decoupled free field theories (free massive 
real fermions) that describe both singlet and triplet sec- 
tors. In contrast in the SU(4) model an analogous de- 
composition does not hold. Indeed, the leading part of 
the bosonized Hamiltonian does not split into two parts 
that account for the six-dimensional (antisymmetric) and 
tenth-dimensional (symmetric) SU(4) irreducible repre- 
sentation: All degrees of freedom strongly interact. In 
the simplest hypothesis, we expect that no phase transi- 
tion occurs between weak and strong couplings but rather 
a smooth cross-over to the plaquette picture described 
above. 

Conclusion - In summary, coming back to the issue 
raised in the introduction, we now have definite evidence 



that the presence of an SU(4) symmetry can indeed have 
very dramatic consequences for lattices in which plaque- 
ttes can form. In the case of the two-leg ladder we have 
shown that there is a spontaneous plaquette formation 
that leads to a degenerate singlet ground state in an oth- 
erwise gapped spectrum. This leads naturally to the con- 
jecture that, for more general lattices, there will be low 
lying singlets, and that such plaquette coverings provide 
a good variational basis to describe them. Work is in 
progress along these lines. 
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